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FOREWORD 

The  report,  'Continuous  Monitoring  of  Dissolved  Oxygen'  is  one 
of  a  series  of  technical  documents  prepared  for  the  Grand  River 
Basin  Water  Management  Study.  The  project  described  herein  was 
undertaken  through  the  Grand  River  Study  Team  at  the  request  of 
the  Grand  River  Implementation  Committee. 

The  material  contained  in  these  reports  is  primarily  technical 
support  information  and,  in  itself,  does  not  necessarily  con- 
stitute policy  or  management  practices.  Interpretation  and 
evaluation  of  the  data  and  findings,  in  most  cases,  cannot  be 
based  solely  on  this  one  report  but  should  be  analyzed  in  light 
of  other  reports  produced  within  the  comprehensive  framework  of 
the  overall  study.  Questions  with  respect  to  the  contents  of 
this  report  should  be  directed  to  the  Co-ordinator  of  the  Grand 
River  Study,  c/o  J.G.  Ralston,  Water  Resources  Branch,  Ministry 
of  the  Environment,  135  St.  Clair  Avenue  West,  Toronto. 
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ABSTRACT 

At  the  time  of  initiation  of  the  Grand  River  Basin  Water 
Management  Study  in  1975,  a  network  of  seven  continuous  moni- 
toring stations  was  established  to  monitor  important  water 
quality  parameters  in  the  central  megalopolis  area  of  the  basin 
at  the  following  locations:  on  the  Grand  River  at  Bridgeport, 
Woolners  Flats,  Blair  and  Glen  Morris,  and  on  the  Speed  River 
at  Canadian  Gypsum  near  Guelph,  Glen  Christie  and  Preston. 
Dissolved  oxygen  and  temperature  are  monitored  continuously  at 
five  stations  using  EIL  (Electronic  Instruments  Ltd.)  instru- 
ment systems.  At  two  other  locations,  NERA  (New  England 
Research  Associates  Inc.)  instrument  systems  record  dissolved 
oxygen,  temperature,  conductivity,  pH  and  redox  potential  at 
half-hourly  intervals. 

The  procedures  involved  in  the  processing  and  storage  of  the 
water  quality  records  of  these  two  instrument  systems  are  out- 
lined in  this  report.  Results  of  statistical  and  probabilistic 
analysis  of  the  dissolved  oxygen  data  collected  during  the 
period  May  1975  -  December  1978  are  also  presented.  The 
results  indicate  that  minimum  daily  dissolved  oxygen  levels 
were  generally  lower  than  the  criteria  levels  of  4  to  5  mg/L 
during  the  months  of  June  and  July  at  all  stations,  with  the 
exception  of  the  Speed  River  at  Glen  Christie  where  the  minimum 
daily  concentrations  were  less  than  4  mg/L  frequently  during 
the  period  May  -  November. 

The  percentage  of  days  in  each  month  in  which  the  minimum 
recorded  dissolved  oxygen  levels  were  in  violation  was  computed 
for  each  station  in  two  ways: 

(1)  as  per  the  former  provincial  criteria  for  warm  water 
biota  (in  which  minimum  concentrations  are  specified); 
and 

(ii)  according  to  the  present  provincial  objectives  for 
warm  water  biota  wherein  the  DO  saturation  levels  (and 
concentrations)  are  specified  as  a  function  of  tem- 
perature. 


vn 


A  comparison  of  the  above  analyses  indicates  that  violations  as 
per  the  present  provincial  objectives  are  generally  less  fre- 
quent than  those  of  the  former  criteria,  except  in  the  case  of 
the  Speed  River  at  Glen  Christie  where  the  opposite  is  true. 
It  is  also  evident  from  this  analysis  that  longer  periods  of  DO 
criteria  violation  occur  under  the  former  criteria,  and  that 
some  periods  of  record  are  in  violation  of  the  old  criteria, 
but  not  the  new  objectives. 
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1.       INTRODUCTION 


1.1      GENERAL 


At  the  time  of  initiation  of  the  Grand  River  Basin  Water 
Management  Study  in  1975,  a  network  of  seven  continuous  moni- 
toring stations  was  established  at  strategic  locations  on  the 
Grand  and  Speed  rivers  in  the  central  megalopolis  area  o*  the 
river  basin.  These  were  instrumented  to  monitor  several  para- 
meters considered  necessary  for  the  diagnosis  of  water  quality 
problems.  Two  more  stations  were  established  during  1973  and 
1979  in  the  lower  megalopolis  area  of  Brantford,  thus  increas- 
ing the  total  number  of  continuous  monitoring  stations  to  nine. 

Technical  Report  No. 11  describes:  (i)  the  selection,  instal- 
lation and  maintenance  of  the  monitoring  systems;  (ii)  pro- 
cessing and  storage  of  data;  and  (iii)  statistical  and  proba- 
bilistic analysis  of  the  data  collected  during  the  period 
1975-1978  at  the  seven  stations  established  during  1975. 

1.2      PURPOSES  OF  CONTINUOUS  MONITORING  OF  WATER  QUALITY 

A  major  objective  of  the  collection  of  continuous  data  on  dis- 
solved oxygen  (DO),  temperature  and  other  water  quality  para- 
meters in  the  Grand  River  basin  is  to  provide  data  required  in 
the  development  and  verification  of  a  continuous  simulation 
water  quality  model.  This  model  when  complete,  will  be  used  to 
evaluate  various  water  management  alternatives  in  the  mega- 
lopolis area  of  the  Grand  River  basin.  Thus,  the  continuous 
monitoring  network  data  will  aid  in  the  long  range  planning  and 
management  of  water  quality  in  the  basin. 

The  continuous  monitoring  data  are  also  useful  in  determining 
the  frequency  of  violation  of  the  provincial  water  quality  cri- 
teria during  selected  seasons   (e.g.,  monthly  intervals). 


These  results  should  be  of  particular  interest  to  those  in- 
volved in  planning  and  surveillance  of  water  quality,  as  well 
as  to  the  general  public. 

The  data  can  be  used  to  identify  water  quality  impacts  attri- 
butable to  an  event-oriented  source  such  as  pollution  due  to 
storm  water  runoff,  as  well  as  to  detect  long-term  trends  in 
water  quality. 

A  potential  application  of  continuous  monitoring  is  the  incor- 
poration into  a  real-time  control  system,  for  management  of 
water  quality  through  proper  scheduling  of  reservoir  releases, 
water  takings  for  water  supply  and  control  of  treatment  plant 
operations.  The  application  of  continuous  water  quality  moni- 
toring in  river  management  is  presented  by  Birch,  et  a^  (1)  and 
Davies  (2). 

At  the  time  of  the  writing  of  this  report,  data  analysis  has 
been  oriented  to  identifying  water  quality  criteria  violations 
and  to  supporting  work  on  the  Grand  River  continuous  simulation 
water  quality  model . 

1.3     SELECTION  OF  SITES 

The  monitoring  network  consists  of  nine  stations  shown  on  the 
basin  map  in  Figure  1.  The  stations  were  selected  based  on  the 
following  rationale: 

Stations  at  the  Canadian  Gypsum  Plant  near  Guelph  on  the  Speed 
River,  and  at  the  Bridgeport  Road  near  Waterloo  on  the  Grand 
River  were  required  to  record  water  quality  at  the  upstream 
boundaries  of  the  megalopolis  area.  The  DO  levels  recorded  at 
these  sites  generally  reflect  the  effects  of  rural  runoff  and 
wastewater  effluents  in  the  upper  reaches  of  the  Grand  and 
Speed  rivers. 
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Figure  1  -  CONTINUOUS  MONITORING  NETWORK 
GRAND  RIVER  BASIN 


The  siting  of  the  remaining  seven  stations,  located  downstream 
of  major  effluent  discharge  points,  was  based  on  estimates  of 
the  downstream  distance  required  for  the  effluent  to  have  mixed 
completely  in  the  entire  cross  section  of  the  river.  The  dis- 
tance required  for  complete  mixing  of  effluent  and  streamwater 
is  usually  known  as  the  "mixing  zone  length".  The  mixing  zone 
lengths  were  estimated  from  the  geometry,  streamflow  and  other 
pertinent  hydraulic  characteristics  of  the  stream  channels. 
Because  of  the  shallow  nature  of  the  streams,  the  distance  re- 
quired for  complete  transverse  mixing  generally  assures  com- 
plete vertical  mixing.  (As  a  rule-of -thumb,  the  distance  re- 
quired for  complete  mixing  in  the  vertical  direction  is  about 
100  times  the  depth;  whereas,  complete  lateral  mixing  is  likely 
to  occur  in  a  distance  of  about  100  times  the  width).  When  the 
mixing  zone  lengths  were  established,  final  site  selection  was 
based  on  practical  considerations  such  as  public  right-of-way 
and  uniformity  of  cross-section  in  a  reasonably  straight  reach. 

The  type  of  instrument  system  installed  at  each  location  is 
shown  in  Table  1.  The  relocation  of  Station  E2,  indicated  in 
Table  1  was  necessitated  by  operational  problems.  A  brief 
description  of  the  systems  is  presented  in  the  next  chapter. 

1.4      PERIODS  OF  RECORD  AND  DATA  ANALYSIS 

Since  the  establishment  of  the  continuous  monitoring  network, 
field  installations  and  operational  procedures  have  improved  to 
the  extent  that  the  sites  can  now  be  kept  in  operation  on  a 
year-round  basis,  except  during  extreme  spring  runoff  condi- 
tions. Figure  2  shows  the  extent  of  the  data  recorded  by 
station  and  month  since  the  installation,  and  indicates  that 
portion  of  the  data  which  has  been  analyzed  and  presented  in 
this  report.  As  stated  before,  the  analysis  presented  in  this 
report  is  confined  to  the  data  collected  during  1975-1978,  at 
the  seven  stations  operational  since  1975. 


TABLE  1.   DETAILS  OF  CONTINUOUS  MONITORING  NETWORK 


Ul 


Location  of 

Station 


Station 
Number 


Grand  River  at 
Bridgeport 

E2 

Grand  River  at 
Woolner  Flats 

E6 

Grand  River  at 
Blair 

N134 

Speed  River  at 
Canadi an  Gypsum 
Plant 

E3 

Speed  River  at 
Glen  Christie 

N135 

Speed  River  at 
Preston 

E5 

Grand  River  at 
Glen  Morris 

E4 

Grand  River 
Upstream  of 

Wilkes  Dam 

E8 

Instrument 
System 


Distance 

above  mouth 
of  Grand  River 
(km) 


Remarks 


Grand  River  at 
Newport  Bridge 


E7 


EIL* 

EIL 

NERA** 

EIL 

NERA 
EIL 
EIL 
EIL 

EIL 


176.8 

166.4 
150.0 
166.0 

163.5 
149.4 
132.8 
108.8 

78.4 


Installed  in  May,  1975  on  the  east  bank  about 

0.5  km  upstream  of  bridge;  relocated  in  April  1977 

on  the  west  bank  about  1  km  upstream  of  the  bridge. 

Installed  in  July  1975 

Installed  in  September  1975 

Installed  in  June  1975  with  housing  on  the  Canadian 
Gypsum  property;  housing  relocated  on  the  Guelph 
WPCP  property  during  Nov.  1978. 

Installed  in  September  1975 
Installed  in  May  1975 
Installed  in  June  1975 
Installed  in  May  1979 

Installed  in  July  1978 


*  EIL  records  DO  and  temperature 

**  NERA  records  DO,  temperature,  pH,  conductivity  and  oxidation  -  reduction  potential 
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O  DATA  COLLECTED   -  NOT  ANALYZED.     a  PART  MONTH  DATA,  ANALYZED.  ■  >20  DAYS  DATA,  ANALYSED. 


Figure  2  -  PERIODS  OP  RECORD  AND  ANALYSIS  OF  DATA 


2.       DESCRIPTION  OF  MONITORING  SYSTEMS 

2.1      INSTRUMENT  SYSTEMS 

Two  types  of  instrument  systems  were  installed  at  the  nine  sta- 
tions of  the  continuous  monitoring  network  in  the  Grand  River 
Basin.  At  seven  locations,  the  monitoring  system  consists  of 
an  EIL  model  15A  meter  (manufactured  by  Electronic  Instruments 
Limited)  in  combination  with  a  Rustrak  model-388  recorder, 
which  records  dissolved  oxygen  (DO)  and  water  temperature  con- 
tinuously on  a  double-strip  chart.  At  two  locations,  NERA 
model-4  monitors  (manufactured  by  the  New  England  Research 
Associates  Inc.),  were  installed.  The  NERA  system  records  dis- 
solved oxygen  (DO),  water  temperature,  conductivity,  pH  and 
oxidation-reduction  redox  potential  on  a  cassette  magnetic 
recording  tape  at  half  hourly  intervals.  Details  of  the  EIL 
and  NERA  instrument  systems  including  description,  operation 
and  maintenance  aspects,  are  provided  in  the  instrument  manuals 
(3,4);  these  are  briefly  reviewed  in  the  following  section. 
Photographs  of  the  instruments  are  shown  in  Figures  3a  and  3b. 

2.2.     DESCRIPTION  OF  EQUIPMENT 

2.2.1    EIL  System 

The  EIL  model  15A  meter  has  two  linear  oxygen  scales  of  0-100% 
and  0-200%,  where  100%  represents  the  concentration  of  oxygen 
in  water  saturated  with  air,  at  a  pressure  of  760  millimeters 
of  mercury.  A  third  scale  is  calibrated  in  temperature  from  -5 
to  +30  C  and  is  used  to  indicate  the  water  temperature*.  The 
dissolved  oxygen  electrode  (model  1511)  gives  an  output  current 
in  proportion  to  the  concentration  of  oxygen  present  in  water, 


*  as  measured  by  the  temperature  electrode  (Model  1512)  in 
the  water. 


Figure  3a  -  EIL  Meter  and  Rustrak  Recorder 
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Figure  3b  -  NERA  Model  4 
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and  this  is  indicated  on  the  ELL  model  15A  meter,  on  a  scale 
directly  calibrated  in  oxygen  percentage  saturation.  Because 
the  output  current  of  the  electrode  varies  markedly  with  temp- 
erature, automatic  temperature  compensation  for  this  variation 
has  been  provided  within  the  temperature  probe  by  the  use  of  a 
thermistor.  The  thermistor  has  a  suitable  temperature  coef- 
ficient which  is  used  to  vary  the  sensitivity  of  the  oxygen 
amplifier.  The  Rustrak  recorder  prints  the  currents  received 
from  the  EIL  meter  through  the  impinging  action  of  its  two 
styli  which  are  driven  by  a  striking  bar  against  the  pressure- 
sensitive  paper.  The  recorder  output  is  a  series  of  dots 
appearing  as  a  continuous  line. 

2.2.2    NERA  System 

The  NERA  model -4  system  consists   of  two  classes  of  equipment: 

(i)  basic  field  equipment,  which  is  taken  to  the  specific 
site  from  which  data  are  to  be  collected;   and 

(ii)  office  equipment,  used  in  the  office  environment  for 
processing  of  the  field  data. 

The  field  equipment  is  portable  and  environmentally  packaged. 
It  consists  of  the  following  subsystems: 

(i)  A  sonde  assembly  which  houses  the  in-situ  type  sensors 
and  protects  them  from  debris  without  interfering  with 
the  environment  to  be  sensed. 

(ii)  A  signal  cable  which  carries  reference  voltages  to  the 
sensors  and  power  to  the  DO  stirrer  motor  and  pream- 
plifiers in  the  sonde  assembly.  The  cable  also 
carries  output  signal  voltages  from  the  sonde  to  the 
signal  conditioner. 


(1ii)  The  signal  conditioner  provides  reference  voltage 
signals  to  those  probes  that  require  them,  and  main- 
tains power  distribution  and  electrical  isolation 
between  probe  sensing  subsystems.  The  signal  condi- 
tioner also  receives  analog  output  signals  from  the 
sonde  assembly,  performs  temperature  compensation 
corrections  to  sensor  signals  as  required  and 
standardizes  the  signals  prior  to  input  into  the 
control  digital  data  logger.  A  multi -scale  meter  on 
the  face  of  the  signal  conditioner  provides  a  quick- 
look  display  of   incoming  data. 

The  control-digital  logger  (C/DDL)  controls  the  timing 
sequence  of  the  overall  system,  determining  when  power 
is  applied  to  probes  that  require  some  settling  time, 
initiating  the  DO  stirrer  sufficiently  ahead  of  data 
scans,  and  initiating  data  scanning  intervals.  The 
C/DDL  performs  analog  to  digital  conversion  of  the 
analog  data  signals  received  from  the  signal  condi- 
tioner. The  digitized  data  signals  are  recorded  on  a 
magnetic  tape  in  such  a  time  sequence  manner  as  to 
permit  subsequent  automatic  handling  and  processing  of 
the  data  by  playback  systems. 

The  present  office  equipment  consists  of  a  DTR-5 
playback/interface  unit  and  software  package  described 
in  Chapter  3.  The  DTR-5  playback  interface  unit 
provides  for  playback  via  a  cassette  recording/  play- 
back device  similar  to  that  used  in  the  C/DDL.  When 
interconnection  to  a  computer  terminal  is  provided, 
the  raw  data  on  the  cassette  tape  is  transmitted  via 
conventional  telephone  lines  to  disk  or  tape  com- 
ponents of  the  Ministry's  computer  facility.  The 
software  provides  for  storage  of  the  raw  data, 
conversion  of  the  raw  data  into  engineering  units, 
error  checks,  storage  of  the  processed  data  and 
display  of  the  processed  data  in  various  formats  as 
described  in  Chapters   3  and  4. 
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Figure  4a  -  EIL  Shelter 


Figure  4b  -  NERA  Shelter 
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2.3  INSTALLATION  OF  INSTRUMENT  SYSTEMS 

Each  EIL  instrument  system  is  housed  in  a  standard  "dog-house" 
which  is  bolted  on  top  of  a  plywood  box  base  (79-cm  x  79-cm  x 
91-cm),  anchor-bolted  to  a  10-cm  thick  concrete  pad  as  shown  in 
Figure  4a.  A  5-cm  diameter  plastic  tubing  is  buried  from  in- 
side the  concrete  pad  to  the  edge  of  the  water,  where  the  tube 
is  connected  to  a  5-cm  cast  iron  conduit.  The  temperature  and 
dissolved  oxygen  probes  are  pulled  through  this  plastic  tube 
and  are  connected  to  the  EIL  meter  inside  the  "dog-house".  In 
the  stream,  the  probes  are  suspended  in  water  with  a  plastic 
boat  float,  and  are  protected  from  algal  fouling  by  an  upstream 
screen  held  vertically  against  the  flow  by  electrical  fence 
rods.  Details  of  an  EIL  monitoring  system  installation  are 
shown  in  Figure  5. 

The  NERA  instrument  system  is  housed  in  a  193-cm  x  173-cm  ARMCO 
shelter  on  a  46-cm  thick  concrete  base  as  shown  in  Figure  4b. 
A  5-cm  polyethylene  conduit  runs  from  the  shelter  to  a  stilling 
well  which  houses  the  sonde  (probe)  assembly.  The  stilling- 
well-sonde  housing  is  constructed  of  a  steel-ridged  culvert 
with  an  attached  15-cm  diameter  plastic  water  pipe.  A  sub- 
mersible pump  near  the  base  of  the  culvert  draws  water  into  the 
culvert  and  discharges  it  downstream  through  a  3-cm  ABS  plastic 
pipe.  Figure  6  shows  the  installation  features  of  the  NERA 
monitoring  system. 

2.4  OPERATION  AND  MAINTENANCE  ASPECTS 

Each  instrument  site  is  visited  once  or  twice  a  week  when  the 
following  checks   and  adjustments   are  carried  out,  as  required: 

(1)  Independent  measurement  of  each  water  quality 
parameter  and  comparison  with  the  meter  reading 
including  resetting  of   reading,    if   required. 
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(2)  Cleaning    slime    and    algal     growth    from    probe    membrane 
and  algal   screen. 

(3)  Recharging  or  replacement  of  battery. 

Both  the  meter  readings  and  independently  measured  values  are 
entered  in  a  table  along  with  the  date  and  time;  this  infor- 
mation is  used  to  correct  the  data  for  calibration  errors 
during  the  processing  of  the  records  as  described  in  Chapter 
3.  Other  checks  and  adjustments  are  also  noted  in  the  table. 
Annual  reports  describe,  in  more  detail,  the  various  aspects  of 
operation  and  maintenance  of  the  continuous  monitors  in  the 
Grand  River  Basin  (5,6). 
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3.      DATA  PROCESSING  PROCEDURES 

3.1  EIL  DATA  PROCESSING 

DO  and  temperature  values  output  by  the  EIL  instruments  are 
recorded  on  a  dual-strip  chart  by  a  RUSTRAK  recorder.  Each 
chart  can  contain  a  maximum  of  about  63  days'  data.  The  DO  and 
temperature  values  recorded  continuously  on  the  strip  chart  are 
digitized  at  predetermined  time  intervals  (i.e.,  converted  from 
chart  co-ordinates  to  numerical  form  on  IBM  punched  cards). 
Briefly,  the  hardware  involved  in  the  digitization  process 
consists  of  an  electrically  charged  (digigrade)  table  top  and  a 
moveable  cursor  which  is  moved  along  a  chart  tracing  which 
translates  chart  points  into  x  and  y  co-ordinates  relative  to  a 
fixed  origin  on  the  table.  The  system  is  linked  to  an  IBM  key- 
punch machine  which  produces  punched  data  from  the  chart  co- 
ordinates. 

The  time  (abscissa)  intervals  are  fixed  for  a  given  chart  and 
are  taken  as  one  hour  for  the  period  of  May  to  November,  when 
diurnal  variations  in  DO  are  significant.  For  the  remainder  of 
the  recorded  year,  an  interval  of  two  hours  is  found  to  be  ade- 
quate as  the  diurnal  DO  variation  is  much  smaller  in  magni- 
tude. These  time  intervals  are  based  on  an  analysis  of  the 
measured  diurnal  DO  profiles,  and  an  optimal  choice  of  interval 
so  as  to  minimize  the  possibility  of  misrepresenting  the  daily 
minimum  and  maximum  in  the  process  of  reduction  of  data  from  a 
continuous   record. 

The  digitized  data  of  the  DO  and  temperature  curves  are  con- 
verted from  RUSTRAK  chart  scale  units  to  real  units  (DO  in  % 
saturation  and  temperature  in  C)  by  the  computer  program  HG 
DOEIL  developed  by  one  of  the  authors  (T.P.H.  Gowda  of  the 
Water  Modelling  Section)  with  assistance  from  J.R.  Eddie  of  the 
Hydrology  and  Monitoring  Section.  The  procedures  involved  in 
the  preparation  of  the  data  for  input  to  the  computer  program 
Are  documented  in  a  manual  (7),  prepared  by  the  Hydrology  and 
Monitoring  Section. 
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Figure  7  shows  a  schematic  sequence  of  the  EIL  data  pro- 
cessing. The  steps  involved  in  the  data  processing  procedure 
are  summarized  below: 

1)  Digitization  of  the  points  on  the  DO  and  temperature 
curves  at  a  predetermined  time  interval  (also  known  as 
sampling  interval). 

2)  Translation  of  the  digital  data  to  DO  and  temperature 
values,   and  time-sequencing  of  each  data  value. 

3)  Application  of  correction  factors  to  each  value  of  DO 
and  temperature  to  account  for  any  errors  recorded 
during  field  inspection  of  the   instruments. 

Two  types  of  errors  may  occur: 

(i)     calibration  differences,   and 

(ii)   improper  resetting  of  the  recorder  needles. 

In  the  former  case,  a  correction  is  applied  on  the  assumption 
that  the  variation  in  the  magnitude  of  error  in  DO  or  tempera- 
ture is  linear  with  time;   the  correction  is  calculated  from 

Hi  -  St?  (ti-tl)  ;  (ti<  ti  <t2) 
t2-  ti 


where  St     ■  estimated  correction  for  parameter  at 


i 


time  t.; 


S      *  difference   between   recorded   parameter 

7 

value  and  independent  measurement  at  end 

of  recording  interval; 
t,     *  beginning  of  recording  period; 
t0  -     end  of  recording  period. 

In  the  case  of  improper  resetting  of  a  recorder  needle,  a 
uniform  correction  is  applied.  A  further  discussion  of  the 
calibration  correction  procedure  is  given  in  Section  3.3. 
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Figure  7  -  eil  data  processing  sequence 
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4)  Storage  of  the  corrected  data  on  a  computer  tape. 

Only  the  initial  hour,  day  and  year  (i.e.,  starting 
time  and  date)  are  stored  on  the  tape  along  with  the 
sampling  interval,  since  the  time  interval  is  constant 
for  any  one  EIL  data  set.  The  meter  calibration  data 
are  also  stored  on  the  tape. 

5)  If  there  are  any  breaks  in  the  chart  recordings,  nega- 
tive numbers  are  substituted  for  data  as  follows: 

-888.00  for  DO;  and 
-999.00  for  Temperature 

Although  this  procedure  can  be  used  for  any  length  of 
time,  it  is  usually  limited  to  90  hours  or  less.  If 
the  break  extends  for  longer  periods,  then  the  data 
are  treated  as  a  separate  set  with  a  new  starting  date. 

6.  When  a  reasonable  number  of  tapes  have  been  created 
for  each  EIL  station  (usually  7  or  more),  the  data  are 
recopied  to  a  single  tape  and  the  older  tapes  erased. 

The  data  stored  on  tape  are  retrieved  by  a  program  in 
which  the  data  points  are  automatically  time-tagged 
using  the  stored  starting  time/date  and  sampling 
interval.  A  sample  output  of  the  DOEIL  program  is 
presented  in  Table  2. 
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TABLE    2:       SAMPLE    OUTPUT    OF    IIGDOEIL    PROGRAM 
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3.2     NERA  DATA  PROCESSING 

3.2.1    Description  of  Computer  Programs 

A  linked  series  of  FORTRAN  programs  are  used  to  process  the 
data  from  NERA  cassette  tapes,  and  to  printout  and  store  the 
results.  These  programs  are  modified  versions  of  programs  ori- 
ginally written  for  use  with  NERA  lake  water  quality  monitors 
maintained  by  the  Lake  Systems  Unit  of  the  Water  Resources 
Branch.  A  complete  documentation  of  the  programs  is  available 
from  the  River  Systems  Unit  (10). 

The  sequence  involved  in  processing  the  NERA  data  records  is 
shown  schematically  in  Figure  8.  The  main  functions  of  the 
individual  computer  programs  are  summarized  below: 

PROGRAM  FUNCTION 

NERA  1  -  selection  of  valid  number  of  channels  (parameters) 

-  selection  of  scale  for  each   channel 

-  read  in  data  from  cassette  tape  to  computer  disk 
storage 

NERA  2  -  read  the  number  of  lines  of  data  to  be  deleted  at 
the  beginning  and  end  of  cassette  tape  record  and 
replace  with  "888888" 

-  check  for  agreement  between  actual  and  expected 
amount  of  data  and  allow  user  cancellation  of  run  if 
discrepancy  is  too  large 

-  selection  of  valid  channels  (parameters),  when  one 
or  more  probes   is  not  functional 

-  select  value  limit  for  each  parameter  for  editing 
out-of-range  values 

NERA  3     -     edit  for  erroneous  characters   in  cassette  tape  record 

-  edit  for  out-of-range  single  values  and  replace  with 
interpolated  values 
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Lgure  8  -  NERA  DATA  PROCESSING  SEQUENCE 
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NERA  4  -  input  field  sheet  data  regarding  monitor  adjustments 
(calibration)  and  breaks  in  record 

-  translate   each   valid   parameter   from   arbitrary 
numbers  to  real  units 

-  apply  calibration  correction  to  data 

-  apply  time/date  tags  to  each  parameter  set 

-  store  edited  data  in  matrix  and  tabular  form  on 
computer  tape 

-  print  out  edited  tabulated  data  (See  Table  3  for 
sample  output. ) 

The  execution  of  these  programs  is  done  through  on-line  (TSO) 
processing  and  requires  about  20  minutes  of  connect  time  per 
cassette  record,  plus  about  one  hour  to  transmit  data  from 
cassette  to  computer  disk  storage.  The  data  on  cassette  tape 
is  read  into  computer  storage  through  a  NERA  Playback/Interface 
unit  model  DTR-5.  This  is  a  110  baud  playback  unit  connected 
between  a  standard  printing  terminal  and  a  telephone  modem, 
which  provides  the  computer  interface.  The  tape  can  also  be 
replayed  off-line. 

3.2.2    Time  Sequencing 

Each  cassette  data  record  starts  and  ends  with  a  row  of  zeros 
("O's")  entered  by  the  field  personnel  changing  the  tape.  When 
the  tape  is  submitted  for  computer  processing,  the  magnetically 
encoded  data  are  transmitted  from  the  cassette  to  a  computer 
disk  file.  The  delimiting  "0"  lines  serve  as  one  check  that 
the  data  set  is  transmitted  in  its  entirety.  The  processing 
routine  also  then  counts  the  number  of  data  in  the  disk  file 
and  prints  out  a  comparison  with  the  number  of  data  expected, 
based  on  the  sampling  interval,  and  the  starting  and  ending 
dates  and  times. 

In  a  normal  processing  run,  the  data  are  then  time-tagged 
sequentially  beginning  at  the  start  date  and  time,  by  adding 
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SAMPLE  OUTPUT  OF  NERA4  PROGRAM 


the  sampling  interval  to  the  time  tag  of  each  subsequent  row  of 
data,  and  the  final  row  of  data  will  then  be  tagged  with  the 
end  date  and  time.  In  the  event  of  a  malfunction  of  the  NERA 
data  logger  during  a  recording  period,  the  field  personnel 
restart  the  cassette  tape  and  punch  in  a  row  of  "O's",  recor- 
ding the  reset  time  and  date  on  the  field  sheet.  During  the 
computer  processing,  this  row  of  "O's"  is  searched  for  by  the 
program  and  when  encountered,  the  subsequent  data  are  re-tagged 
starting  with  the  reset  time. 

3.3  COMMENTS  ON  AUTOMATED  DATA  PROCESSING 

The  present  system  employed  in  the  processing  of  continuous 
monitor  data  results  from  a  good  co-operative  effort  between 
the  field  and  office  staff,  and  has  produced  a  ^ery  reliable 
data  base  for  the  period  1975-1978.  However,  some  minor  prob- 
lems occasionally  arise  which  require  special  handling  of  the 
data: 

(a)     Data  Calibration  Correction 

As  discussed  above,  when  instrument  decalibration  occurs,  a 
linear  correction  over  time  can  be  applied.  The  possibility 
exists  however,  that  the  deviation  is  non-linear,  either 
because  the  deviation  occurs  intermittently,  or  because  the 
deviation  varies  with  the  time  of  day.  In  the  case  of  an 
intermittent,  or  sudden  change  in  DO,  for  example  that  result- 
ing from  a  sudden  fouling  of  probes  by  a  drifting  mat  of  algae, 
an  EIL  chart  record  will  clearly  show  the  effect;  the  personnel 
digitizing  the  data  can  then  interpolate  values  on  the  chart 
during  the  problem  period.  When  this  is  done,  the  calibration 
procedure  is  bypassed  in  the  data  processing.  When  a  similar 
event  occurs  during  the  recording  of  NERA  data,  the  absence  of 
a  chart  record  makes  such  sudden  anomalous  DO  changes  less 
obvious.  A  computer  routine  to  detect  such  periods  is  being 
considered. 
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Every  effort  is  made  to  check  the  accuracy  of  the  data  from 
both  types  of  monitors,  and  although  automated  data  processing 
is  necessary  to  handle  the  large  amounts  of  data  collected,  the 
system  is  flexible  enough  that  data  are  corrected  manually  when 
necessary. 

(b)     NERA  Data  Timing 

EIL  data  are  recorded  originally  continuously  on  a  chart  and 
can  be  visually  checked  for  time  continuity.  However,  NERA 
data  are  not  time-tagged  until  the  cassette  data  record  is 
processed.  This  has  led  to  problems  on  a  few  occasions. 

Over  a  period  of  months,  the  sampling  interval  in  one  NERA 
monitor  gradually  lengthened  from  30  to  34  minutes.  This  would 
have  led  to  large  discrepancies  in  time-tagging  of  data  if  it 
had  not  been  noticed  during  the  data  processing.  This  occur- 
rence, together  with  the  occasional  loss  of  some  data  during 
transmission  of  the  cassette  data  to  computer  disk,  led  to 
modifications  to  the  program  NERA  2  which  now  ensures  that  the 
entire  set  of  data  is  received  on  disk,  and  is  properly  time- 
tagged. 

Secondly,  data  logging  equipment  malfunctions,  or  running  out 
of  tape,  has  caused  breaks  in  the  data  record.  The  program 
NERA  4  has  the  capability  of  re- tagging  data  points  according 
to  information  supplied  on  the  field  sheets,  in  this  event- 
uality. However,  when  both  of  the  above  mentioned  problems 
also  occur  in  the  same  recording  period  (which  occurred  on  one 
occasion),  it  is  difficult  to  detect  the  problem  and  produce 
reliably  timed  data. 

The  following  improvements  in  NERA  data  processing  procedures 
are  therefore  suggested: 

(i)  Specification  and  implementation  of  a  continuous  means  of 
verifying  the  time  accuracy  of  the  NERA  data  logging 


-  26  - 


system.  Any  sporadic  malfunctions  or  gradual  change  in 
the  scanning  interval  would  thus  be  immediately  obvious  to 
field  personnel  and  could  be  rectified  without  any  delay. 
The  present  system  of  having  the  amount  of  data  verified 
by  the  data  processing  system  ensures  the  timing  accuracy 
of  the  data,  but  there  is  an  unavoidable  time  lag  in 
reporting  any  discrepancies  to  the  field  personnel. 

(ii)  Further  data  processing  program  modifications  to  ensure 
that  NERA  data  are  corrected  appropriately  where  cali- 
bration corrections  must  be  made. 

3.4      ACCURACY  OF  PROCESSED  DATA 

The  accuracy  of  the  processed  data  (i.e.,  outputs  of  the 
HGDOEIL  and  the  NERA  computer  programs)  is  dependent  on  the 
magnitudes  of  the  potential  instrumental  errors  as  well  as 
manual  errors.  A  brief  description  of  the  results  of  limited 
investigations  on  these  error  sources  is  presented  below. 

(a)     Instrumental  Errors: 


The  specifications  for  the  probes  used  in  the  EIL  meter  state 
that  the  accuracy  of  temperature  measurements  is  within 

+  1.5  C,  and  that  the  DO  readings  are  accurate  to  within  +  5% 

o 
of  the  scale  readings  for  solutions  within  10  C  of  a  pre- 
viously selected  temperature  (3).  The  specifications  for  the 
temper-  ature  probes  of  the  NERA  instrument  system  state  that 
the  over-  all  accuracy  of  the  output  data  signal  is  +  0.2  C 
for  the  temperature  range  -5  C  to  25  C,  +  0.4  C  for  the 
range  25  C  to  45  C,  and  the  overall  accuracy  of  the  meter 
reading  is  stated  to  be  the  same  as  that  of  the  output  data 
signal  +  0.25°C  (4).  The  output  data  signal  of  the  NERA  DO 
probe  is  specified  to  have  an  overall  accuracy  of  +  2%  of  the 
reading;  the  overall  accuracy  of  the  meter  reading  is  specified 
to  be  {+  2%  of  the  reading  +  0.5%  of  the  range). 
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Investigations  were  carried  out  to  determine,  under  field  con- 
ditions, the  differences  between  the  EIL  DO  probe  output  and 
the  RUSTRAK  recorder  chart  readings.  The  absolute  differences 
of  the  temperature  outputs  were  found  to  be  within  0.5°C  at 
all  seven  EIL  stations,  except  for  differences  of  up  to  3  C 
at  times  during  early  spring  and  late  fall.  The  DO  readings 
were  found  to  be  accurate  to  within  +  1%  saturation. 


Results  of  investigations  at  the  two  NERA  sites  showed  that  the 
temperature  recordings  were  accurate  to  within  +  0.15  C, 
whereas  the  accuracy  of  the  DO  readings  was  within  +0.1  mg/L. 

(b)  Manual  Errors 

These  errors  apply  to  the  EIL  data  and  are  generated  during  the 
digitization  of  the  chart  records.  Errors  are  introduced  when 
the  cursor  is  not  centered  exactly  at  the  desired  points  on  the 
DO  or  temperature  curve.  In  order  to  determine  the  magnitude 
of  this  error,  separate  samples  (each  consisting  of  24  data 
points)  were  selected.  A  comparison  of  the  DO  values  of  the 
processed  data  (i.e.,  computer  outputs)  with  the  corresponding 
chart  readings  showed  an  absolute  mean  difference  of  about  \h%, 
with  a  range  of  -2%  to  +3%;  Generally,  a  faster  rate  of  digi- 
tization was  associated  with  a  larger  discrepancy.  Under 
normal  conditions,  an  accuracy  of  1  to  2%   (  +  )  is  achieved. 

(c)  Total  Error 

Based  on  the  above  facts,  the  processed  data  is  likely  to  be 
accurate  to  within  +  1  C  for  temperature  and  +  4%  saturation 
for  dissolved  oxygen. 
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4.       STATISTICAL  AND  PROBABILISTIC  ANALYSIS  OF  DATA 

4.1      GENERAL  DETAILS  OF  DATA  ANALYSIS 

Since  the  inception  of  the  continuous  monitoring  networks  in 
the  Grand  River  Basin  in  1975,  a  vast  amount  of  data  has  been 
generated  on  DO,  temperature,  oxidation-reduction  potential,  pH 
and  conductivity.  As  an  example,  with  a  sampling  interval  of 
one  hour  a  record  of  one  year  would  result  in  the  generation  of 
8760  values  of  each  parameter  at  each  monitoring  station.  In 
order  to  summarize  the  data  in  a  compact  form  without  a  signi- 
ficant loss  of  the  most  important  characteristics  of  the  data 
records,  the  statistical  and  probabilistic  data  analysis  proce- 
dures, described  in  the  following  subsections,  have  been  used. 
This  approach  is  also  consistent  with  the  requirements  of  the 
continuous  monitoring  data  for  calibration  and  verification  of 
the  Grand  River  water  quality  simulation  model.  As  stated 
earlier,  this  analysis  includes  the  data  generated  at  the  seven 
stations  installed  in  1975  for  the  recording  period  May  1975  - 
December  1978. 

A  schematic  summary  of  the  statistical  and  probabilistic  data 
analyses  is  presented  in  Figure  9.  Brief  descriptions  of  the 
data  analysis  programs  are  presented  below. 

4.1.1    Daily  Statistics 

By  considering  a  daily  period  of  24  hours  from  the  midnight  of 
one  day  to  the  midnight  of  the  next  day  (EST),  the  daily  sta- 
tistics were  computed  by  the  program  DAILY,  which  retrieves  an 
EIL  or  NERA  dataset  and  provides  the  following  output  summary: 

(i)   daily  mean  value  of  temperature  (°C)  and  DO  (mg/L); 
(ii)  daily  minimum  value  of  temperature  (°C)  and  DO  (mg/L), 


and  times  of  occurrence; 
daily  maximum  value  of  t< 
and  times  of  occurrence. 


(iii)  daily  maximum  value  of  temperature  (  C)  and  DO  (mg/L), 
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Figure  9  -  CONTINUOUS  MONITORS  RAW  DATA  ANALYSES 


A  sample  output  of  the  DAILY  program  is  presented  in  Table  4. 
The  statistical  summary  tables  for  the  period  May  1975  - 
December  1978,  are  on  file  and  are  available  on  request. 
Although  the  program  is  set  up  to  print  output  summaries  of 
only  temperature  and  DO  at  present,  it  has  the  capability  of 
producing  similar  output  for  the  additional  parameters  of  con- 
ductivity, pH  and  redox  potential,  recorded  by  the  NERA  moni- 
tors. 

4.1.2    Probabilistic  Analysis 

Part  of  the  output  from  the  DAILY  program  was  used  to  derive 
exceedance  probability  plots  for  DO.  The  percentage  proba- 
bility of  exceedance  of  each  value  was  calculated  by  ranking 
the  data  for  a  given  month  in  descending  order  of  DO  concen- 
tration, and  using  the  equation 


where  p  =  percentage  probability  of  exceedance  of  a  given 
value; 
m  =  ranking  of  each  data  value;  and 
N  -  total  number  of  data  points  in  the  month. 

The  probabilistic  analysis  was  performed  with  minimum  daily  and 
mean  daily  DO  concentrations  for  all  seven  stations  for  the 
recording  period  May  1975-December  1978.  The  computations  were 
carried  out  by  using  the  program  SORT. 

A  preliminary  evaluation  of  the  arithmetic  and  logarithmic 
probability  plots  of  the  minimum  and  mean  daily  DO  data  indi- 
cated a  straight  line  relationship  on  the  latter  plot;  this 
suggests  that,  in  general  the  data  conform  to  a  log-normal 
distribution.  Based  on  this,  the  results  were  plotted  on  loga- 
rithmic probability  paper,  and  straight  lines  of  best  fit 
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were  drawn  through  the  data  points.  In  some  cases,  the  data 
were  found  to  deviate  from  the  straight  line  of  best  fit  in  the 
lower  and  higher  range  of  probabilities  (generally  less  than 
20%  and  greater  than  90%).  This  is  thought  to  result  from  a 
lack  of  homogeneity  of  the  DO  data  at  its  more  extreme  values. 
However,  these  plots  were  found  to  be  useful  for  characteriza- 
tion of  the  DO  for  comparison  between  stations  and  identifica- 
tion of  seasonal  trends,  as  well  as  giving  a  clear  indication 
of  the  station  DO  characteristics  at  a  glance. 
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4.2      GRAND  RIVER  AT  BRIDGEPORT  -  STATION  E2 

The  probability  distribution  plots  of  mean  daily  DO  and  minimum 
daily  DO  for  each  month  at  station  E2,  the  upstream  control 
station  on  the  main  Grand  River,  are  presented  in  Figure  10. 
It  should  be  noted  that  data  were  not  available  for  the  period, 
December  1976-May  1977.  Generally,  the  difference  between  the 
minimum  and  mean  daily  values  during  the  winter  months  is  very 
small  compared  with  the  difference  during  the  rest  of  the 
months.  The  DO  levels  were  mainly  above  5  mg/L,  with  some 
exceptions  during  the  months  of  June  and  July  when  the  DO 
levels  were  less  than  4  mg/L  for  up  to  20%  of  the  time.  The 
break  in  the  probability  plot  for  July  1976  resulted  from 
changes  arising  from  a  very  large  increase  in  the  streamflow 
rate. 


34 


BRIDGEPORT  1975 


tr> 


OS 


td 


U 


Legend: A  -  Mean  Daily  D.O.;  o  -  Minimum  Daily  D.O. 


u 


DEC.  7  5 


i  lnnimil  i  ■  t  il  i  Mil  nnlmijtiiil  n  1 1  hllill  1 1  ll  1 1 1  I  I  I  I    llllll 

10    20   30   40  50   60  70   30   ^0  93 


Probability  of  Occurrence  -  % 

FIGURE  10  -  PROBABILITY  DISTRIBUTION  OF  D.O. 
LEVELS  AT  STATION  E2  -  BRIDGEPORT 


-35- 


Legend : 


BRIDGEPORT  1976 
Mean  D^iiy  n.o.j  °  -  Minimum  naiiy  n.n 


X 


14 


10 


MAY   7fi 


10 


ei 


JUNE  76 


10     20   30  40   50  60  70    80    90 
Probability  of  Occurrence,  % 


99 


FIGURE  10  -  CON'T. 


36- 


BRIDGEPORT  1976 
Legend:  A  -  Mean  Daily  D.O.;  o  -  Minimum  Daily  D.O. 


V 


2 
O 


E- 
2 
W 

u 

2 
O 

o 


S   10 


10 
2 


NOV,  7G 


I  UllllllUllllfllllllllllllllllllll|llll|||||||||||i|[|[|||| 

10     20   30  40   50  60  70    80    90 


20   30  40   50  60  70    80    90 
To]  ability  of  Occurrence,  % 


LJUl^ 


FIGURE  10  -  CON'T. 


BRIDGEPORT  1977 


Legend: A  -  Mean  Daily  D.O.;o  -  Minimum  Daily  D.O. 


10 


2    5 


JUNE  7' 


M 


10 


AUG.  77 


w 


i  '  1 1  ■  ill imiiifl  i  ■  I il i  n il ii 


DEC.    77 

ilmilniiliiiiliinniiiliiil  i  i  i    inn 


10  20       30       43    50    60       70       80  90 


9  9 


Probability   of   Occurrence,    % 
FIGURE    10-    CON'T 


-38- 


BRIDGEPORT  1978 


Legend:^  -  Mean  Daily  D.O.  ;  o  -  Minimum  Daily  D.O. 


53 


m 


o 


15 

10 

_        A 

A  A 

~$~~ 

A    4 

r""^" 

^> -* 

©_, 

5 

/ 
/ 

6 

MAY    78 

u 


O    5 


13 

10 

5 

ft 

-fi — V          f                                                JUIiY  78 

iTiiiTitiiiiiiiliiiiliiiiliiiiliiiilitiiliiiiliiiiiiiiiliiiiii  i  inn 

10    20    30  40   50  60  70   80     90 
Probability  of  Occurrence,  % 


99 


FIGURE  10  -  CON'T. 


-39- 


BRIDGEPORT  1978 


Legend:  A  -  Mean  Daily  D.O.;  q-   Minimum  Daily  D.O. 


31 

a 


10 


10 


0 


AUG.  7J 


—  1  iiiiihiiiiiiiriiiiiliiiiliiiiliiiiliiiiliiiiTiiiiiimliiiii  1  1  1111 


10     20   30  40   50  60  70   80     90 
Probability  of  Occurrence,  % 


99 
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4.3      GRAND  RIVER  AT  WOOLNER  FLATS  -  STATION  E6 

Figure  11  shows  the  monthly  probability  distributions  of  mini- 
mum and  mean  daily  DO  at  station  E6  located  on  the  main  branch 
of  the  Grand  River  between  Waterloo  and  Kitchener.  No  data 
were  available  for  the  period  May-June  1977.  The  data  points 
are  seen  to  follow  the  log-normal  distribution  in  a  majority  of 
the  cases;  however,  in  some  of  the  plots,  e.g.,  Oct.  77,  there 
is  a  scatter  of  the  data  points  in  the  lower  probability 
range.  DO  levels  were  generally  high,  with  the  exception  of 
July  1977  and  June-September  1978. 
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WOOLNER  FLATS  1978 


Legend:  A  -  Mean  Daily  D.O. ;  o  -  Minimum  Daily  D.O. 


en 


m 


w 


u 


0 

■n «T~ 

$    ®~ 

"W 

-§    °^ 

6 

OCT. 

m 

14 
z 

o       10 

o 


14 


10 


llllllllllllllllllll 


10 


20   30  40   50  60  70   1ST)    9CT 
Probability  of  Occurrence  -  I 


FIGURE  11  -  CON'T. 


-47- 


4.4     GRAND  RIVER  AT  BLAIR  -  STATION  N134 

Data  generated  between  November  1977  and  October  1978  at  the 
NERA  monitoring  station  N134,  located  between  Kitchener  and  the 
Grand-Speed  confluence,  have  been  included  in  the  probability 
plots  of  mean  and  minimum  daily  DO  levels,  shown  in  Figure  12. 
From  these  plots,  the  DO  levels  are  seen  to  be  greater  than  4 
mg/L  for  most  of  the  recording  period.  Only  July  and  August 
show  any  periods  of  DO  less  than  5.0  mg/L. 
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FIGURE  12  -  CON'T. 
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FIGURE  12  -  CON'T. 
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4.5      SPEED  RIVER  AT  CANADIAN  GYPSUM  PLANT  NEAR 
GUELPH  -  STATION  E3 

The  DO  probability  plots  for  station  E3,  the  upstream  station 
on  the  Speed  River,  are  presented  in  Figure  13.  The  data  show 
that  the  DO  levels  were  over  4  mg/L  during  the  entire  period  of 
record;  the  plots  show,  however,  that  up  to  40%  of  the  days  had 
minimum  DO  levels  below  5  mg/L  in  July  and  August  of  1975,  1977 
and  1978.  Although  this  station  is  located  at  the  downstream 
boundary  of  the  City  of  Guelph,  the  DO  levels  are  not  likely  to 
have  been  adversely  affected  by  urban  storm  water  runoff  due  to 
the  very  short  travel  times  involved. 
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FIGURE  13  -  PROBABILITY  DISTRIBUTION  OF  D.O. 
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4.6      SPEED  RIVER  AT  GLEN  CHRISTIE  -  STATION  N135 

The  DO  probability  plots  at  the  NERA  station  N135,  located  be- 
low the  Guelph  water  pollution  control  plant  (WPCP),  are  shown 
in  Figure  14  for  the  period  June  1976-December  1978,  with  the 
exception  of  February  and  March  1977.  The  minimum  DO  plots 
contrast  with  those  of  other  stations  in  that  the  lines  of  best 
fit  are  generally  steeper  than  those  for  the  mean  DO  plots. 
This  corresponds  to  a  greater  variation  of  within-day  DO  range 
than  for  other  stations.  Generally,  during  the  period  May- 
November  of  each  year,  the  DO  levels  are  seen  to  be  very  low, 
being  frequently  lower  than  4  mg/L.  The  daily  mean  DO  plots 
show  less  than  5  mg/L  for  up  to  60%  of  the  time  during  several 
months. 

The  break  in  the  plot  of  September  1976  data  is  probably  due  to 
a  washout  of  algal  biomass  resulting  from  a  large  increase  in 
flow  rate  in  mid-month.  Plots  for  several  other  months  illu- 
strate similar  deviations  from  the  linear  form  indicating  a 
sensitivity  of  DO  at  this  station  to  physical  changes  in  the 
river. 
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4.7     SPEED  RIVER  AT  PRESTON  -  STATION  E5 

Figure  15  shows  the  probability  distribution  plots  of  the  mean 
and  minimum  DO  levels  at  station  E5,  located  near  the  mouth  of 
Speed  River,  for  the  period  May  1975  -  December  1978.  Gener- 
ally, the  data  follow  the  log-normal  distribution.  From  the 
plots,  it  can  be  seen  that  the  DO  levels  were  very  low  during 
May- September,  minimum  daily  DO  being  lower  than  4  mg/L  on  up 
to  60%  of  the  days  in  some  months.  Data  was  available  for  June 
through  August  1978,  but  has  been  omitted  from  the  plots  be- 
cause of  a  probe  malfunction  during  this  period. 
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4.8     GRAND  RIVER  AT  GLEN  MORRIS  -  STATION  E4 

Station  E4  is  located  at  the  downstream  edge  of  the  central 
Megalopolis  area  of  the  river  basin.  Probability  distribution 
plots  of  mean  and  minimum  DO  levels  for  the  period  May  1976  - 
December  1978  are  presented  in  Figure  16.  There  were  no  occur- 
rences of  minimum  daily  DO  less  than  4.0  mg/L  in  May  of  1975, 
1976  and  1978;  however,  there  were  several  occurrences  of  mini- 
mum daily  DO  less  than  4.0  mg/L  in  June,  July  and  August  of 
1976  and  1977.  DO  levels  of  between  4.0  and  5.0  mg/L  are 
comnon  during  this  summer  period. 
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4.9      PROBABILITY  OF  OCCURRENCE  OF  DAILY  MINIMUM  DO  LEVELS 

The  DO  probability  plots  shown  in  figures  10  to  16  indicate 
that  the  minimum  daily  DO  levels  were  below  5  mg/L  during  the 
months  of  June  and  July  at  several  stations.  From  these  proba- 
bility plots,  station  summary  plots  were  derived;  these  plots, 
presented  in  Figure  17,  show  the  percentages  of  days  in  each 
month  during  which  the  minimum  DO  levels  did  not  exceed 

4  mg/L,  and  5  mg/L,  at  each  station. 

In  the  main  branch  of  the  Grand  River  between  Bridgeport  and 
the  confluence,  the  minimum  daily  DO  levels  were  occasionally 
lower  than  4  mg/L  during  the  months  of  May,  June  and  July, 
whereas  the  DO  levels  were  lower  than  5  mg/L  on  up  to  80%  of 
the  days  during  May-September.  The  DO  levels  required  for  the 
protection  of  warm  water  biota  are  only  marginally  being  met 
during  the  months  of  May-September. 

In  the  Speed  River  between  Guelph  and  the  confluence,  the  mini- 
mum DO  levels  are  lower  than  4  mg/L  more  frequently.  At  Glen 
Christie  (Station  N135),  the  minimum  daily  DO  levels  often 
approach  zero,  and  are  observed  to  be  lower  than  4  mg/L  for  up 
to  100%  of  the  time.  From  the  Glen  Morris  plot,  it  can  be  seen 
that  the  minimum  daily  DO  levels  were  less  than  or  equal  to 

5  mg/L  on  up  to  60%  of  the  days  during  May-September. 
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4.10     VIOLATION  ANALYSIS 

Two  types  of  violation  analysis  were  carried  out  in  order  to 
interpret  the  results  according  to  the  traditional,  as  well  as 
the  present  provincial  DO  criteria.  Brief  descriptions  of  the 
violation  analysis  procedures  are  presented.  A  comparison  of 
the  results  of  these  two  types  of  violation  analysis  is  dealt 
with  in  subsection  4.10.3. 

4.10.1   Violation  of  Traditional  DO  Criteria 

Formerly,  the  dissolved  oxygen  criteria  were  stated  in  terms  of 
minimum  concentrations  of  DO,  irrespective  of  water  tempera- 
ture. The  provincial  DO  criteria  for  warm  water  biota  were 
stated  as  follows: 

The  dissolved  oxygen  (DO)  concentration  should  be  above  5 
mg/L  at  all  times,  except  that  in  certain  situations  con- 
centrations may  range  between  5  and  4  mg/L  for  short  inter- 
vals within  any  24-hour  period  provided  that  water  quality 
is  favourable  in  all  other  respects.  (8) 

In  order  to  interpret  the  DO  records  using  these  criteria, 
daily  violation  summaries  of  DO  concentrations  were  obtained  by 
using  the  computer  program  VIOLAT.  This  program  prints  out  a 
daily  summary  of  the  number  of  hours  during  which  the  DO  level 
was  at  or  below  a  specified  concentration  level,  viz:  1,2,  ...  9, 
and  10  mg/L  or  greater.  The  same  tabulation  is  also  shown  as  a 
percentage  of  time  during  which  the  DO  levels  were  at  or  below 
the  specified  concentration  levels.  A  typical  output  of  the 
VIOLAT  program  is  shown  in  Table  5. 
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4.10.2   Violation  of  Present  Provincial  DO jriteria 

The  present  provincial  water  quality  objectives  are  stated  as  a 
function  of  temperature  and  are  specified  for  two  types  of 
biota:  warm  water  biota  and  cold  water  biota.  It  was  decided 
to  apply  the  objectives  for  the  warm  water  biota  in  this  analy- 
sis. The  DO  criteria  for  warm  water  biota  are  as  follows: 

At  no  time  should  dissolved  oxygen  concentrations  be  less  than 
the  values  specified  below: 


Temperature,  C 

0 

5 

10 

15 

20 

25 

Dissolved 

Oxygen 

Concentration 

% 

saturation 

47 

47 

47 

47 

47 

48 

— j 

mg/L 

7 

6 

5 

5 

4 

4 

In  order  to  obtain  DO  violation  summaries  as  per  these  cri- 
teria, a  computer  program  VIONEW  was  developed.  This  program 
prints  out  the  percentage  time  of  violation  of  the  DO  levels 
for  the  above  range  of  temperatures  for  each  day,  as  well  as 
the  total  time  of  violation  during  that  day.  A  typical  output 
of  the  VIONEW  program  is  presented  in  Table  6.  In  this  analy- 
sis it  is  assumed  that  the  precise  criterion  is  the  percent 
saturation  value  and  not  the  concentration  value  (also  shown  in 
the  objectives),  since  the  latter  is  an  approximation  of  the 
former. 

The  tabulated  violation  summaries  for  the  period  May  1975  - 
December  1978,  obtained  from  the  VIOLAT  and  VIONEW  programs, 
are  on  file  and  are  available  on  request. 
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TABLE  6:   SAMPLE  OUTPUT  OF  VIONEW  PROGRAM 


4.10.3   Comparison  of  Violation  Analysis  Results 

In  order  to  compare  the  violations  of  the  foregoing  types  of  DO 
criteria,  average  percentage  times  of  DO  violations  in  each 
month  were  calculated  as  follows: 

(1)  From  the  VIOLAT  program  output,  the  daily  average  per- 
centage of  data  violating  the  4  mg/L  criterion  was 
computed  for  each  station. 

Similarly,  the  daily  average  percentage  of  data  vio- 
lating the  5  mg/L  criterion  was  computed  for  each 
stati  on. 

(3)  The  daily  average  percent  of  data  violating  the  pre- 
sent DO  objectives  for  warm  water  biota  in  each  month 
was  computed  from  the  VIONEW  program  output. 

These  results  are  compared  in  Figure  18.  In  some  months,  the 
former  DO  criteria  would  be  in  violation,  whereas  the  new 
criteria  are  not  being  violated.  Generally,  violation  of  the 
new  criteria  is  more  frequent  than  violation  of  the  4  mg/L 
level,  but  is  less  frequent  than  violation  of  the  5  mg/L  level. 

At  all  stations  except  the  Speed  River  station  at  Glen  Christie 
(N135),  the  percentages  of  violations  of  a  5  mg/L  criterion 
would  be  higher  in  comparison  to  violations  of  the  new  provin- 
cial objectives.  In  addition,  months  occur  where  violations 
are  in  evidence  for  the  former  criterion  of  5  mg/L,  but  not  for 
the  new  objectives  (e.g.,  summer  1978  at  station  E2-Bridge- 
port).  However,  at  the  Speed  River  NERA  Station  (N135),  the 
percentages  of  violations  of  the  new  provincial  objectives  are 
higher  than  those  of  the  former  criteria;  the  months  in  which 
violations  occur  are  almost  the  same  in  both  cases. 
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FIGURE  18  -  AVERAGE  TIME  OF  VIOLATION  OF  DO  CRITERIA 
AT  VARIOUS  STATIONS  IN  THE  GRAND  R.  AND 
SPEED  RIVER   DURING  1975-1978 
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FIGURE    18    -(continued) 


5.      SUMMARY  AND  CONCLUSIONS 

Technical  Report  No.  11  deals  with  the  continuous  monitoring  of 
water  quality  at  seven  stations,,  operational  since  1975,  in  the 
central  megalopolis  area  of  the  Grand  River  basin.  At  five 
stations,  DO  and  temperature  are  monitored  continuously  by  EIL 
instrument  systems;  at  two  other  locations,  DO,  temperature, 
conductivity,  pH  and  oxidation-reduction  potential  are  recorded 
at  half  hourly  intervals  by  NERA  instrument  systems.  The 
results  of  the  statistical  and  probabilistic  analyses  of  the  DO 
data  collected  at  the  seven  stations  during  the  period  1975-78, 
as  presented,  appear  suitable  for  calibration  and  verification 
of  the  continuous  simulation  water  quality  model,  which  is 
being  developed  for  application  to  the  Grand  River  system. 

An  analysis  of  the  magnitude  and  duration  of  violations  of 
provincial  warm  water  DO  criteria  by  recorded  DO  levels  has 
been  carried  out.  The  violation  analysis  is  based  on  former 
provincial  DO  criteria  (in  which  minimum  concentrations  are 
specified),  as  well  as  the  present  provincial  water  quality 
objectives  wherein  the  DO  requirements  for  aquatic  biota  are 
specified  as  a  function  of  temperature.  The  DO  data  record  was 
evaluated  in  terms  of  these  two  types  of  criteria  and  a  compar- 
ison made. 

The  following  conclusions  summarize  the  results  of  the  analyses 
presented  in  the  report: 

1.  The  DO  levels  in  the  Grand  River  at  Bridgeport  (Station  E2) 
were  generally  above  5  mg/L  with  some  exceptions  during  the 
months  of  June  and  July  when  the  minimum  DO  levels  were 
below  4  mg/L  on  up  to  20%  of  the  days. 

2.  At  Woolner  Flats  in  the  Grand  River  (Station  E6),  the  DO 
concentrations  were  high  relative  to  provincial  objectives, 
with  some  exceptions  during  July  and  August  1977,  and 
June-September  1978. 
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3.  In  the  Grand  River  at  Blair  (Station  N134),  the  DO  levels 
were  generally  higher  than  4  mg/L  during  most  of  the 
recording  period. 

4.  With  the  exception  of  July  1975,  the  DO  levels  in  the  Speed 
River  at  Guelph  (Station  E3)  were  greater  than  4  mg/L; 
however,  for  some  summer  months,  the  daily  minimum  DO 
concentrations  were  lower  than  5  mg/L  on  over  30%  of  the 
days.  This  station  had  the  most  acceptable  DO  conditions, 
relative  to  provincial  objectives,  of  any  of  the  stations 
with  a  four  year  data  record. 

5.  In  the  Speed  River  at  Glen  Christie  (Station  N135),  the 
minimum  DO  levels  were  frequently  lower  than  4  mg/L  during 
the  period  May-November.  In  1978,  the  daily  minimum  DO 
values  were  below  5  mg/L  for  the  period  May  to  October,  and 
in  July  and  August  the  minimum  values  were  below  4  mg/L. 
This  station  had  the  most  unacceptable  DO  conditions,  rela- 
tive to  present  provincial  objectives,  of  any  of  the 
continuous  DO  monitoring  stations. 

6.  At  the  mouth  of  the  Speed  River,  at  Preston  (Station  E5), 
the  daily  minimum  DO  levels  were  lower  than  4  mg/L  on  up  to 
60%  of  the  days  during  some  months. 

7.  On  the  Grand  River  at  Glen  Morris  (Station  E4),  the  daily 
minimum  DO  concentrations  were  less  than  4  mg/L  during 
May- July  on  up  to  30%  of  the  days,  and  less  than  5  mg/L 
during  the  same  months  on  up  to  60%  of  the  days. 

8.  A  comparison  of  the  DO  violations  of  former  provincial 
criteria  against  present  provincial  DO  objectives  for  warm 
water  biota  indicates  that  more  violations  of  the  fixed  5 
mg/L  criterion  occur  than  violations  of  the  new  objec- 
tives. This  is  true  for  all  stations  except  station  N135 
on  the  Speed  River  at  Glen  Christie,  where  the  opposite  is 
true.  Most  stations  therefore  show  months  during  which 
violations  occur  of  the  former  criterion,  but  not  the 
present  objectives. 
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